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Abstract 

Computations of divergence and vorticity are made by the tri- 
angular method for six 100-nb layers (1C00-400 mb) for those obser- 
vation times (ll) when wind data is available from eight stations 
in the BOMEX region during the period 11-28 July 1969 (Phase IV). 

The results are compared to cloudiness from satellite pictures. 
Contemporary correlation coefficients between vorticity and di- 
vergence in the six layers ranged from -.19 to +.26. The corre- 
lation in the lowest layer (1000-900 mb) (where Ekman theory pre- 
dicts a large (~ -1.0) negative correlation) was -.19. The vorticity 
and divergence values had the same sign in approximately 50/o of the 
sample. Despite the low values for the correlation coefficients, 
the comparison of the results to satellite pictures is encouraging 
in that convergence and cyclonic vorticity are associated with re- 
gions of maximum cloudiness. When algebraic means are computed 
within each of the six layers, the results indicate the vorticity 
has the same order of magnitude from 1000 to 400 mb while the 
divergence decreases by an order of magnitude above 800 mb. The 
divergence and vorticity were of the same order of magnitude from 
1000 to 800 mb. Computation of mean magnitudes showed no such de- 
crease in the order of magnitude of the divergence above 800 mb. 

All triangles were subjectively classified as either disturbed, 
undisturbed or neither. Comparison of the mean values of diver- 
gence and vorticity within layers by the Student T s t-test revealed 
there were significant (l$ level) differences between the vorticity 
results for disturbed and undisturbed conditions and the results 
for divergence indicate significant differences at the 2 £ level in 
the lowest layer and at the 10-20^ levels in the next two layers 
(900-700 mb) thereby supporting the CISK hypothesis. 
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INTRODUCTION 

The Barbados Oceanographic and Meteorological Experiment (BOMEX) 
was carried out from 3 May to 28 July 1969 east of the island of Bar- 
bados (13,1 N, 59.5 W). The entire experiment was divided into four 
phases and it is phase IV (11-28 July 1969) with which this thesis 
is concerned. This phase of BOMEX was designed to explore large 
convective systems in the tropical Atlantic. The overall plans, 
description and some early results of BOMEX activities are described 
in Davidson (1968), Cook (1969), Kuettner and Holland (1969), U.S. 
Dept, of Commerce (1969), Holland (1970), Friedman and McFadden 
(1970) and BOMEX Bulletins numbers 1-9. 

There were five "stationary" ships plus the island of Barbados 
to take observations during phase IV. The geographic arrangement is 
shown in figure 1. In addition, there were numerous aircraft avail- 
able which were equipped to take meteorological observations. During 
the period 11-28 July 1969 there were six synoptic-scale disturbances 
which passed through the BOMEX array as described by Fernandez- 
Partagas and Estoque (1970) , one of which was classified as a tropi- 
cal depression (22-27 July). The remaining disturbances were classi- 
fied as follows: 3 weather systems, 1 complex weather situation and 

1 middle-level weakening cyclonic circulation (Fernandez-Partagas 
and Estoque, 1970). This crude method of classifying tropical dis- 
turbances suggests a definite lack of understanding of the distur- 
bance-producing mechanisms in the tropics. 
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Fig. I. Geographic arrangement of fixed stations during 
BO MEX , phase J3L (Scale , 1/2 " = 1° ) . 
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In looking for causes of tropical disturbances, Charney and 
Eliassen (1964) and Charney (1969) have proposed conditbnal insta- 
bility of the second kind (CISK) as a mechanism for establishment 
of low-level water-vapor convergence and a resulting upward motion 
at the top of the boundary layer, thus leading to a release of latent 
heat above the boundary layer which provides the source of energy 
for maintaining a disturbance. CISK requires an initial increase 
in the surface relative vorticity as the initiating impulse to set 
up the boundary-layer convergence. Gray (1968) found a . .strong 
association of tropical disturbance and storm development with 
synoptic-scale surface relative vorticity . . The idea of CISK 

as a cause of disturbances in the tropics has met with wide approval, 
but observational evidence of such a mechanism is not available; in 
fact, how to obtain observations to confirm the hypothesis is not 
readily apparent to the author. Sargeant and Ruttenberg (1970) have 
included the following as one of the "Specific Atmospheric Objec- 
tives" of an international tropical experiment which will probably 
be conducted in 1974: "In particular, to provide sufficient data 

to verify or reject present formulations of the CISK mechanism, and 
to extend such formulations". 

Manabe et, al. (1970) have attempted to reproduce the general 
circulation features of the tropics from a numerical model for the 
month of January and obtained quite good results. They simulated 
the convection process by means of "moist-convective adjustment", 
a process which has the effect of transferring heat energy from 
lower layers to upper layers, thus creating a cold core in the lower 
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part of the troposphere and a warm core in the upper troposphere. 

This process was initiated in the model whenever a layer became 
saturated and the lapse rate in the layer became super-moist 
adiabatic . 

In an attempt to extend present ideas on tropical convective 
phenomena, it was decided to use the BOMEX data to look at distur- 
bances in some detail. Although the distance between stations (380- 
530 km) precludes investigation of small-scale activity, it does pro- 
vide the necessary data for investigation of the disturbance-scale 
(i.e., sub-synoptic scale) as well as for comparing differences in 
conditions between disturbed and undisturbed situations. The rela- 
tion of the mesoscale structure of the tropics in disturbed vs. un- 
disturbed conditions is now being recognized as one of the most 
crucial problems to solve in the tropics (Zipser, 1970). 

The array of ship and land stations in the BOMEX area is ar- 
ranged conveniently for application of the triangular method of com- 
puting divergence, vertical motion and vorticity as proposed by 
Bellamy (1949). A modification of Bellamy’s method was proposed by 
Graham (1953) but since the data stations in BOMEX are "fixed 11 the 
Bellamy method is appropriate. The results thus obtained are most 
likely to be representative of the centroid of the triangle (Bellamy, 
1949) . Figure 2 shows the breakdown of the BOMEX area into eleven 
triangles with centroids marked by (2) . Two stations have been 
added to expand the area of coverage for which data is available 
twice daily. The major disadvantage of the Bellamy method is the 
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rainier 




Fig. 2. Triangular divisions with expanded BOMEX array. 

Numbers at centroids correspond ta numbers assigned 
to that triangle for identification .( Scale : 1/2" = 1° ) 
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assumption of a linear variation of the wind between vertices of the 
triangle, an assumption which may yield misleading results in cer- 
tain cases as shown by Zipser (1965). Successful application of 
this method has been made by Byers and Rodebush (1948) and Byers 
and Hull (1949) in thunderstorm studies. Lateef (1967) computed 
divergence, vertical motion and vorticity fields over a grid cover- 
ing the Caribbean for a 3-day period during which a low-level 
easterly wave passed through the grid. He used the pressure- 
differentiated form of the continuity equation to compute vertical 
velocities to offset the difficulty of obtaining a large non-zero 
value for the net divergence in a column. Various other authors 
who have made computations of divergence, vertical motion and vor- 
ticity in the tropics utilizing the continuity equation include 
Landers (1955) , Arnason (1955) , Endlich and Mancuso (1963) , Kyle 
(1970) and Arnason et. al. (1963). The computations of Landers 
(1955) utilized the modified Bellamy method proposed by Graham 
(1953) while the method used by Arnason et. al. (1963) consisted 
of extracting wind components at gridpoints from streamline and 
isotach analyses. Kyle (1970) computed vertical motion using a 
latitude-longitude grid of monthly mean values of the wind compo- 
nents. Endlich and Mancuso (1963) performed their calculations by 

the triangular method on triangles varying in area from 22000 to 
2 

920000 km in the Caribbean area. In contrast, the triangles shown 

2 

in figure 2 range in size from approximately 71000 to 105000 km . 
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APPLICATION 

In order to be able to synoptically analyze the results, only 
those times when winds for all eight stations were available were 
used. This limited the data to 00Z and 12Z observations because 
that is the only data available for Raizet and Chaguaramas which 
is found in the Northern Hemisphere Data Tabulations (NHDT) . This 
means a possibility of 18 days at twice per day which would yield 
36 data sets. However, so much of the ship wind observations 
were missing that the end result was 11 data sets, and they aren’t 
all quite complete. The observation times finally selected for 
this study are: 14/OOZ, 14/12Z, 18/00Z, 20/00Z, 21/12Z, 22/12Z, 
25/12Z, 26/OOZ, 26/12Z, 28/00Z and 28/12Z. Tables 3 to 5 in Appen- 
dix B summarize the identification of the triangles and the loca- 
tions of the centroids. 

In order to interpret the results in comparison with disturbed 
or undisturbed situations, a method of classifying the triangles 
as disturbed or undisturbed is needed. Such a method is available 
by means of enlargements of ATS-III and ESSA-9 satellite pictures. 

By tracing the cloud cover in the BOMEX region onto acetate and then 
projecting the result with an overhead projector, enlargements of 
the satellite picture can be obtained to almost any scale. For this 
study the pictures were enlarged to a scale of 1" = 1° of latitude 
(results presented later are at a scale of 1/2 M = 1° of latitude). 
With a copy of figure 2 at the same scale superimposed on the copy 
of the satellite picture, a subjective classification of U (undis- 
turbed) , D (disturbed) or N (neither disturbed nor undisturbed) was 
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made for each triangle at each observation time depending on the areal 
extent of cloud cover for each triangle. Allowances were made for 
cases when the satellite picture was not at the same time as the ob- 
servations. The criteria used for the classifications were as 
follows: U, less than 10% of area covered by clouds; N, less than 
50% of area covered by clouds but more than 10%; D, 50% or more of 
area covered by clouds. Out of a total of 121 such classifications 
(11 triangles for each of the 11 observation times used) there were 
52 = U, 27 = N and 42 = D. The results of this classification are 
listed in Appendix B, table 6. 

Appendix A describes the method of computing divergence and 
vorticity. In addition to the triangle parameters, the data needed 
for the computations are the wind direction and wind speed profiles 
for each of the stations. For this study the winds were obtained 
at 50-mb intervals from 1000 to 400 mb and then vectorially averaged 
by 100-mb layers according to the trapezoidal rule, thus resulting 
in six layers, each containing a vector-mean wind. The vector- 
mean winds were then used to compute the mean divergence and mean 
vorticity for each of the six layers of each of the eleven triangles. 
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DATA UTILIZED 

The data utilized for the computations were the u and v wind 
components at 50-mb intervals from 1000 to 400 mb. The data for 
Barbados, Raizet and Chaguaramas were used exactly as given in 
the NHDT. Data for the five "stationary" ships were complicated 
by the fact that the ships were actually not stationary due to the 
failure of the mooring lines for all five ships before the begin- 
ning of phase IV (BOMAP, 1971). Since it was planned to use the 
processed rawinsonde observations from the n A 0 " (BOMEX Bulletin 
No, 9, 1971) data reduction process performed at NASA's Mississippi 
Test Facility, this presents the problem of correcting the winds 
for ship's motion, A search through ship's logs and course and 
speed tables provided by BOMAP (BOMAP, 1971) revealed that there 
were several inconsistencies in regards to the necessary correction. 
(For example, ship's course before, during and after could not be 
verified when compared to gyro headings found in the BOOM data 
[BOMAP, 1971]), Fortunately, however, it appears that this cor- 
rection is very small (because tabulated ship speeds are small) for 
all ships except ROCKAWAY. The tabulated ship speeds for ROCKAWAY 
were as large as 10.5 knots. Due to the uncertainty in such a 
large correction it was decided to use the data sent by ROCKAWAY 
via teletype with the idea that the observers on ROCKAWAY made die 
necessary corrections to the data before sending it. 

An additional problem with ROCKAWAY was the inability to ac- 
quire the balloon for 6 to 10 minutes after launch (BOMAP, 1971), 
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as evidenced by the missing winds in the lowest 5000 feet or so. 

To alleviate this problem the teletype winds were broken down into 
u, v components and plotted on a vertical logarithmic scale using a 
standard tropical atmosphere (Schacht, 1946; and Riehl, 1954). The 
winds at the lowest reported level were then linearly connected to 
the surface wind and the necessary winds were read off at 50-mb 
intervals . 

Except for specific cases to be noted later, the "A 0 " data 
was used for the remaining ships (at present there is no choice 
but to use it for Mount Mitchell). The u, v winds were taken di- 
rectly from the M A 0 n data at 50-mb intervals. In a few cases 
winds were missing at some of the 50-mb intervals (in particular, 
for DISCOVERER) and in such cases the needed data was obtained by 
applying techniques similar to that described above for ROCKAWAY. 

The only corrections made to this data were for RAINIER which had 
a voltage reference problem. Comment cards acquired from BOMAP 
on the performance of each individual sounding indicate that the 
computed wind directions for RAINIER are in error by 140-220°. This 
is particularly evident if wind directions computed from !, Ao M data 
are compared to wind directions reported in the teletype data. Con- 
sequently the error reported by BOMAP was applied to all the wind 
directions for RAINIER as computed from the "A 0 " data. Even 
though no correction was tabulated for RAINIER at 26/12Z, a cor- 
rection of 140° was applied based on soundings at times before and 
after 26/12Z and based on comparisons with the teletype data. 
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Teletype data was also substituted for "A 0 " data in the following 
three cases for the reasons specified: 

21/12Z Oceanographer ; Many winds were missing in the "A 0 " 
data but complete wind data was available in teletype form. 

22/12Z Oceanographer; There were no winds in the f, A© M data 
but complete wind data was available in teletype form. 

22/12Z Discoverer; The "A 0 M winds were missing at many 
levels but complete wind data was available in teletype form. 
There was no data at all for OCEANOGRAPHER at 28/12Z but there was 
"A©" data for 28/10Z and 28/1343Z. In order to complete the data 
set at 28/12Z, the winds at the two off times were averaged and 
used in the computations to obtain an approximation to the values 
for 28/12Z. 

A final comment regarding corrections to the island data is 
necessary. No terrain effects have been considered for Barbados, 
Raizet and Chaguaramas. As pointed out by Lilly and LaSeur (1956) 
the nature of a correction at Chaguaramas is very complicated. 

This was also discussed by Zipser (1965) who made "ad hoc n correc- 
tions to the Chaguaramas winds. Zipser’ s correction to the winds 
at Raizet amounted to increasing the reported mean wind speed 
(surface - 850 mb) by 17%. For this study, no corrections have 
been applied to the reported winds for any of the island stations. 
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RESULTS 

Some of the results of this study are presented in this section 
in the form of synoptically analyzed charts of the BOMEX region. 
Vorticity and divergence for layers 1 and 2 (1000-900 mb and 900- 
800 mb) for five synoptic times are presented in the figures 3-12. 
Results for all synoptic times used and all layers are presented 
in tables 7-12 in Appendix B. In addition, figures 3-12 show de- 
pictions of cloud cover obtained from ATS-III and ESSA-9 satellite 
pictures for comparison to the analyzed results which are presented. 
These depictions were obtained from the satellite pictures as des- 
cribed in the Application section. 

The figures presented in this section cover two specific 
cases: 1) 14 July 1969; Weather system (Fernandez-Partagas and 

Estoque, 1970) which moved through the northern half of the BOMEX 
region from east to west and 2) 25-26 July 1969; Tropical depression 
(Fernandez-Partagas and Estoque, 1970) which emerged from the ITCZ 
and moved northwest through the BOMEX region. 

A comparison of the satellite picture depiction at 14/00Z (keep 
in mind the satellite picture is for 13/1817Z) with the analysis of 
divergence and vorticity shows nothing very striking, perhaps due 
to the relative inactivity in the BOMEX region. Perhaps the most 
noticeable feature is the convergence which seems to be coming into 
the region from the east. At 14/12Z the convergence has extended 
slightly farther into the BOMEX region and appears to be somewhat 
associated with the region of cloudiness which is approaching from 




Fig. 3. Analyzed results of divergence and 
v orticity for I4/00Z , layer I % with 
cloud cover for 13/1817 2* 






Fig. 5: Analyzed results of divergence and vortici>y 
for 14/12 Z , layer I , with cloud cover 
for 14/ 11362. 
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the east. Likewise there is a region of cyclonic vorticity which 
appears to be associated with the cloud system and anticyclonic vor- 
ticity out ahead of it. 

Comparing the 25/12Z results (the satellite picture reproduc- 
tion is for 25/1532Z) shows, the ITCZ to be associated with the region 
of strongest convergence. The stipled regions in the cloud depic- 
tion correspond to regions on the picture which were very bright. 
Comparison with the vorticity analysis shows the ITCZ region to be 
associated with the cyclonic vorticity; however, the orientation of 
the vorticity lines appears to be transverse to the ITCZ. 

For comparison of results at 26/00Z it becomes necessary to use 
a 25/1826Z satellite picture. Unfortunately there isn f t much dif- 
ference between the 25/15Z and 25/18Z pictures but the difference 
in the results for divergence and vorticity is considerable. The 
region of strongest convergence is still associated with the region 
of maximum cloudiness but now there is also a strong (^4 x 10 “*) 
center of cyclonic vorticity which appears to be located in the 
ITCZ region. 

At 26/12Z the pattern has changed considerably. There is a 
broad expanse of cloud cover which covers almost’ the entire BOMEX 
region. There is a strong (-3.5 x 10 “*) center of convergence in 
the center of the region and the small region of divergence (13N, 

61W) agrees well with a clear region on the satellite picture. The 
vorticity results are not as good except for the cyclonic vorticity 
region in the northern half of the BOMEX region. 




Fig. 7: Analyzed results of divergence and vorticity 
for 25/122 , layer I , with cloud cover 
for 25/1532 2- 





Fig. 8: Same as Figure 7 except layer 2 . 





Fig. 9-. Analyzed results of divergence and vorticity 
for 26/00 2 , layer I , with cloud cover 
for 25/1826 2. 
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(or 26/1233?. 
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In order to compare these results with those predicted by Ekman 
theory, a correlation coefficient between divergence and vorticity 
was computed. Correlations were computed for all layers to look 
at differences between layers. To minimize the effects of round- 
off error, all values were multiplied by 10“* before computing any 
means or sums. The results of this are as follows (% with same 
sign means percent of sample where vorticity and divergence had 
same sign) : 



Layer (mb) 


Sample Size 


% Same Sign 


Correlation 
f icient 


1000-900 


117 


43.6 


-.19 


900-800 


117 


55.5 


1 

o 

u> 


800-700 


113 


58.4 


+ .26 


700-600 


113 


56.6 


+.15 


600-500 


109 


49.5 


-.05 


500-400 


96 


50.0 


+.08 



These values in the lower layers are somewhat disappointing from 
what one would expect from Ekman theory. Nevertheless it must be re- 
membered that this is a linear correlation coefficient and as such is 
merely a measure of the scatter about a best-fit line of the values of 
divergence and vorticity. Another possibility is that there is a lot 
of "noise 1 * in the data which was used. It was anticipated that by 
using layer-mean winds in the computations, instead of winds at dis- 
crete levels, much of the unwanted noise or inconsistencies would be 
minimized. Perhaps the mixing of "Ao" data with teletype data cE w?ll cS 
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linearly adjusting the ROCKAWAY’s low-level winds contributed to 
data noise. The low values for the correlation coefficient could 
have been anticipated from the higher-than-expected percentage 
(50%) of cases where the divergence and vorticity had the same sign. 

The possibility of computing divergence and vorticity on one 
scale with data contaminated by local small-scale winds may also be 
a problem. It merely means the assumption of a linear change in the 
wind between vertices of the triangle does not apply. By using layer- 
mean winds, however, the effects of such wind regimes should be 
minimized. 

The profiles of the mean values of divergence and vorticity 
provide an interesting result which is summarized in table 1. (The 
six samples which were used to obtain the results in table 1 are 
contained in Appendix B, tables 7-12). The significant fact which 
is immediately obvious is that the mean (algebraic) divergence is 
an order of magnitude greater below 800 mb than it is above while 
the magnitude of the mean (algebraic) vorticity is about the same 
from 1000 to 400 mb. The sample size is not very large even in 
the lowest layers where there is more data available, but the maxi- 
mum values of cyclonic vorticity and convergence do appear to occur 
in the lowest layers. Furthermore, in the mean, the vorticity and 
divergence are of opposite sign and are of the same order of mag- 
nitude from 1000 to 800 mb. 

Table 1 also contains the values of the mean magnitude of 
divergence and vorticity within each of the six layers. When looked 
at in this manner the magnitude of the divergence does not decrease 
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Table 1. Summary of mean values of divergence and vorticity. 
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by a factor of 10 above 800 mb. Comparison of these results with 
those of Endlich and Mancuso (19 o3j p. 21) show a tendency for the 
magnitude of the vorticity to be 5-20^ larger than theirs up to 
500 mb while the magnitude of the divergence is 20-60fo larger than 
theirs. Above 500 mb the magnitudes of both divergence and vorticity 
given in table 1 are less than those of Endlich and Mancuso (1963). 

Similar computations of algebraic means and mean magnitudes 
were carried out for each of the three categories of degree of 
disturbance (U, N and D) and the results are presented in table 2. 
Although this results in three small samples, particularly for 
category N, the results between categories do appear to differ 
dramatically in several respects. The mean (algebraic) vorticity 
is quite large and cyclonic from 1000 to 600 mb in category D 
whereas it is weak and anticy clonic in category U. The mean 
vorticity is also cyclonic from 1000 to 600 mb in category N but 
not as large as in category D. The variance of the vorticity in 
category D is quite large and increases from 1000 to 700 mb even 
though the algebraic mean remains practically constant. The mean 
(algebraic) divergence also differs considerably between categories 
U and D. In category U there is very slight convergence but it is 
an order of magnitude less than the convergence from 1000 to 800 mb 
in category D. The results for category N appear to be intermediate 
to those of categories U and D (as would be expected) except for the 
convergence in layer 3* 

To determine if the results in table 2 show any significant 
differences between categories, the Students t-test (Panofsky and 



Category U 



Layer 
.(mb) _ 


Sample 

Size 


JL 


Vari- 

ance 


Ifl 


w 


Vari- 

ance 


lv-\v| 


1000-900 


51 


-.05 


1.12 


.89 


-.08 


.93 


.74 


900-800 


51 


-.26 


1.76 


1.11 


0.00 


1.09 


.83 


800-700 


50 


-.13 


1.09 


.83 


-.04 


1.12 


.83 


700-600 


50 


.06 


1.20 


.89 


.01 


1.72 


1.00 


600-500 


49 


-.19 


1.37 


.87 


.08 


2.01 


1.01 


500-400 


43 


-.75 


1.10 


1.01 


-.13 


1.19 


.89 



Category N 



1000-900 


24 


.56 


1.06 


.98 


-.31 


I.67 


1.04 


900-800 


24 


.26 


1.43 


.92 


-.24 


2.92 


1.33 


800-700 


23 


.12 


1.70 


.90 


-.26 


2.39 


1.17 


700-600 


23 


.12 


2.44 


1.01 


.02 


1.53 


.82 


600-500 


23 


-.23 


1.34 


.93 


.38 


1.45 


.92 


500-400 


21 


-.20 


.90 


.84 


.13 


1.01 


.83 



Category D 



1000-900 


42 


.94 


2.44 


1.46 


-.61 


1.43 


1.05 


900-800 


42 


.98 


2.51 


1.5^ 


-.27 


1.93 


1.08 


800-700 


40 


.98 


3.10 


1.66 


.17 


.96 


.78 


700-600 


40 


.66 


2.69 


1.4l 


0.00 


1.10 


.84 


600-500 


37 


-.02 


1.92 


1.13 


-.08 


1.36 


• 93 


500-400 


32 


-.20 


1.01 


.75 


-.04 


.94 


.71 



Table 2. Summary of mean values of divergence and vorticity _5 

within categories of degree of disturbance. Units are 10 . 
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Brier, 1958) was performed on the computed means of vorticity and 
divergence by layers between the three categories. If the results 
are significant they would be expected to occur between categories 
U and D. 

The results of this test were as follows: 







Categories 


U vs. D 




Layer 


if 


$ Chance 
that U=D 


t 7 .„ 


$> Chance 
that U=D 


1 


3*574 


Cl 


2.333 


2 


2 


4.067 


Cl 


1.071 


10-20 


3 


3.662 


cl 


.932 


10-20 


4 


2.073 


<3 


.024 


>45 


5 


.617 


25-30 


.538 


25-30 


6 


2.204 


<3 


.363 


30-40 






Categories 


U vs. N 




Layer 


tl 


$ Chance 
that U=N 




$ Chance 
that U=U 


1 


2.313 


C3 


.837 


20-25 


2 


1.591 


5-10 


.753 


20-25 


3 


.859 


10-20 


.702 


20-25 


4 


.194 


40-4-5 


.023 


>45 


5 


.124 


45 


.884 


10-20 


6 


1.912 


2-5 


.918 


10-20 



Categories N vs. D 



Layer 


ts 


$ Chance 
that N=D 




$ Chance 
that N=D 


1 


1.041 


10-20 


.941 


10-20 


2 


1.926 


2-5 


.074 


>45 


3 


2.005 


1-2 


1.313 


5-10 


4 


1.270 


10-20 


.046 


>45 


5 


.596 


25-30 


1.451 


5-10 


6 


.029 


>45 


.619 


25-30 



The results are evidently very significant. During disturbed 
conditions there is significantly more cyclonic vorticity than 
during undisturbed conditions. Furthermore there is significantly 
more convergence in the lowest layer in disturbed versus undisturbed 
conditions and the results for convergence in the 900 to 700-mb 
layers could easily be called meteorologically significant. The 
results between categories U and N and categories N and D occassion- 
ally show significant differences but this is. probably caused by the 
N-saraple being a mixture of both U and D values as well as being 
a small sample. It is interesting to note that there is apparently 
no significant difference between samples for convergence in the 
700 to 600-mb layer because all three samples were essentially 
non-divergent according to the algebraic mean (see table 2). 

The results between categories U and D are in direct support 
of the CISK hypothesis (Charney and Eliassen, 1964) in that there 
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is increased cyclonic vorticity and increased convergence in the 
lowest layers during disturbed conditions. In addition the results 
presented here indicate the increased cyclonic vorticity is 
present up to 600 mb whereas the increased convergence occurs up 
to 800 mb. 
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CONCLUSIONS 

From the results presented in this paper, it appears that the 
triangular method of computing divergence and vorticity in the 
tropics yields acceptable results. The contemporary correlations 
of vorticity and divergence in the lowest layers were disappointing, 
however the good agreement which was obtained when comparing sate- 
llite pictures to analyzed results was encouraging. Yfhen the results 
were divided into categories depending on whether they occurred dur- 
ing disturbed or undisturbed conditions as determined by cloud cover 
from satellite pictures, there were significant differences between 
disturbed and undisturbed conditions for both vorticity and diver- 
gence. There was a significant increase in cyclonic vorticity and 
convergence which is in support of the CISK hypothesis (Charney and 
Eliassen, 1964). Algebraic mean values of divergence and vorticity 
for the entire sample had the same order of magnitude below 800 mb 
but the divergence decreased by an order of magnitude above 800 mb. 
Furthermore the mean (algebraic) vorticity retained the same order of 
magnitude from 1000 to 400 mb. 

It is suggested that additional computations be made of the 
nature described in the previous sections but without the emphasis 
on those times when all eight stations reported. Rather the empha- 
sis should be shifted so as to obtain a very large sample of data 
to determine if the results in tables 1 and 2 remain valid. A more 
objective method of classifying the results into categories of undis- 
turbed or disturbed is needed. In addition it is not possible from this 
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study to conclude the reason for the divergence and vorticity having 
the same sign 50% of the time in the lowest layers, a result which 
is contradictory to Ekman theory. 

Classification of the triangular regions as disturbed or un- 
disturbed might be improved by combining the method used in this 
study with a method which is more objective. Such a method is avail- 
able using ©e. the equivalent-potential temperature (Rossby, 1932 
and Hess, 1959). Garstang et. al. (1967) found that higher values 
of occurred at all levels when the atmosphere was disturbed. 

In some preliminary reduction of the BOMEX data for 12-15 and 
23-28 July 1969 ( 168 soundings) the author found the vertically- 
averaged value of © e (surface to 400 mb) to vary from 325. OK to 
342. IK, the higher values occurring during disturbed conditions. 

The advantage of using 0 e is that it combines the effects of tem- 
perature and moisture into one parameter. Furthermore such an ob- 
jective method could be used by itself for times when no satellite 
pictures are available. Such an approach might yield results which 
show more significant differences between the mean divergence and/or 
vorticity in disturbed and undisturbed conditions. Significant 
differences (5% level according to Student’s t-test) were not 
found in this study although the differences between the undis- 
turbed versus disturbed divergence in the lowest three layers were 
significant at 15-35% levels. 

Due to the difficulties encountered with the ’’Ao" data format, 
it would seem that the teletype data is still the best (most reliable) 
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data source, even though it T s been almost two years since the ter- 
mination of BOMEX. Since the observers on board the ships are "pro- 
fessionals” at their jobs, it seems reasonable that any corrections 
which need to be made to the data (e.g. ship motion) would be made 
by the observers before sending the results of a launch. It is 
unfortunate that detailed records are not readily available which 
reflect what corrections, if any, were made by the observers be- 
fore sending the teletype message. 

Since no corrections were made to the winds at the island sta- 
tions, it is suggested that a method for adjusting the winds repor- 
ted at these stations be adopted before conducting additional re- 
search with triangles affected by these stations. Corrections were 
suggested by Zipser (1965) for Raizet and Chaguaramas ; however the 
"correction" used by him at Chaguaramas was "ad hoc". 
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APPENDIX A 



The Computation Procedure 

Numerous triangles can be formed within the BOMEX array as 
shown in figure 2 which can be used to calculate horizontal con- 
vergence ( ) as shown by Bellamy (1949): 




where I Vi I is the observed wind speed, hu = the altitude of the 
triangle through vertex i, d_^ = the observed wind direction and 

= the azimuth of the side of the triangle opposite vertex i. 

With M in (m sec and lu in (m) , the units of become sec \ 

Bellamy (1949) describes construction of an overlay to make the 
necessary computations or, as an alternative, to construct tables 
which can be entered with wind direction aid wind speed. However, 
with the advent of the electronic computer direct computations using 
(1) are simplified. 

Computations of vorticity can also be performed using (1) if 
d^ is increased by 90° (Bellamy, 1949), such that 



f - t 

J i - « Wl 



( 2 ) 



Since data were used at constant pressure surfaces, the compu- 
tation of vertical motion with the continuity equation is straight- 
forward. 



c>CO 






( 3 ) 
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Integrating leads to 

CO, = cO, - V-VAP 



/t/VWv 

where the symbol C *) implies a vertical average. 

Using the boundary condition that 

CO =0 at p = 1000 mb 

enables computation of vertical velocities by successive integra- 
tions of (3) . ( A review of various methods of computing vertical 
motion are contained in Panofsky [1946]). An additional assump- 

A'VVW* 

tion that V-\\/ - v-v was also made. 



(4) 



(5) 
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APPENDIX B 

Table 3. Latitude and longitude of data stations: 



Station 


Latitude 


(N) 


Longitude (W) 


i. 


Oceanographer 


7.5 




52.7 


2. 


Mount Mitchell 


10.5 




56.5 


3. 


Discoverer 


13.0 




54.0 


4. 


Rockaway 


15.3 




56.6 


5. 


Rainier 


17.5 




54.0 


6. 


Barbados 


13.1 




59.5 


7. 


Raizet 


16.3 




61.5 


8. 


Chaguaramas 


10.7 




61.6 


Table 4. 


Designation of 


triangles : 






Triangle No, 


Vertices of 


the triangle 




i 


Oceanographer , 


Chaguaramas, Mount Mitchell 




2 


Oceanographer , 


Mount 


Mitchell , Discoverer 



3 

4 

5 

6 

7 

8 
9 

10 

11 



Mount Mitchell, Chaguaramas, Barbados 
Mount Mitchell, Barbados, Discoverer 
Discoverer, Barbados, Rockaway 
Mount Mitchell, Barbados, Rockaway 
Mount Mitchell, Rockaway, Discoverer 
Discoverer, Rockaway, Rainier 
Rockaway, Raizet, Rainier 
Barbados, Raizet, Rockaway 
Chaguaramas, Raizet, Barbados 



44 



APPENDIX B 

Table 5. Latitude and longitude of triangle centroids: 



Triangle No. 


Latitude(N) 


Longitude (W) 


1 


9.5 


56.9 


2 


10.3 


54.4 


3 


11.4 


59.2 


4 


12.2 


56.6 


5 


13.8 


56.7 


6 


13.0 


57.5 


7 


12.9 


55.7 


8 


15.2 


54.9 


9 


16.3 


57.4 


10 


14.9 


59.2 


11 


13.4 


60.9 



Table 6: Subjective evaluation of disturbed nature of each triangle for each data time. 



45 




Triangle No. 14/00Z 14/12Z 18/OOZ 20/00Z 21/12Z 22/12Z 25/12Z 26/OOZ 26/12Z 28/00Z 28/12Z 

1 ’ 1.22 -0.71 0.02 -2.11 -2.36 2.93 0.02 -0.04 -0.57 
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TABLL 8.Results for 900-800 mb. Description the same as for table 



Triangle No. 14/OOZ 14/12Z 18/00Z 20/00Z 21/12Z 22/12Z 25/12Z 26/OOZ 26/12Z 28/00Z 28/12Z 

1 -0.32 0.24 -0.13 0.07 1.96 2.21 -0.70 -0.42 -0.01 
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TABLE 9. Results for 800-700 mb. Description the same as for table 



Triangle No. 14/00Z 14/12Z 18/00Z 20/00Z 21/12Z 22/12Z 25/12Z 26/00Z 26/12Z 28/OOZ 28/12Z 

1 -1.72 0.70 -1.53 -2.25 -0..00 0.99 2.24 1.56 1.64 
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TABLE 10- Results for 700-600 mb. Description the same as for table 



Triangle No. 14/00Z 14/12Z 18/00Z 20/00Z 21/12Z 22/12Z 25/12Z 26/00Z 26/12Z 28/00Z 28/12Z 
1 -1.66 0.34 1.69 . 0.19 0.62 0.63 1.62 0.42 0.87 
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Results for 600-500 mb. Description the same as for table 



Triangle No. 14/00Z 14/12Z 18/OOZ 20/00Z 21/12Z 22/12Z 25/12Z 26/00Z 26/12Z 28/00Z 28/12Z 

1 0.20 1.38 1.16 2.73 1.35 0.52 0.60 -0.22 
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TABLE 12. Results for 500-400 mb. Description the same as for table 
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